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We study the validity of the large-N equivalence between four-dimensional SU(N) lattice gauge 
theory and its momentum quenched version — the Quenched Eguchi-Kawai (QEK) model. We 
have found strong evidence that this equivalence does not hold in the weak-coupling regime (and 
thus in the continuum limit). This is based on weak-coupling analytic arguments and Monto- 
Carlo simulations at intermediate couplings with 20 < N < 200. Since detailed expositions of 
our arguments, methods and results have already appeared in Phys. Rev. D78:034507 (2008) and 
Phys. Rev. D78:074503 (2008), we attempt here to give a more intuitive explanation of our results. 
The breakdown of reduction that we find is due to a dynamically generated correlation between 
different Euclidean components of the gauge fields. 
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1. Introduction 

QCD simplifies in the 't Hooft limit of a large number of colors, and as a result it has been 
a long-standing goal to understand the properties of the theory in that limit [jl]], including on the 
lattice [0]. In this paper we reconsider an alternative to conventional large volume simulations, 
namely the use of large-Af volume reduction to single-site models [||, f|, |5| |6|] (see also the related 
Ref. [0]). This allows one, in principle, to study very large values N ~ 0( 100 — 400) with modest 
resources. In this paper we choose to study one of the variants of the original Eguchi-Kawai (EK) 
volume reduction, namely the quenched Eguchi-Kawai (QEK) model. Our motivation is two- 
fold: First, it is the only single-site model whose equivalence with large-Af QCD has yet to be 
thrown in doubt (in contrast to twisted Eguchi-Kawai model — see Refs. Second, we use it as 
a tool to gain experience with single-site large-Af models, before turning to the theory with adjoint 
fermions [0]. 

Since a detailed discussion of this work, including a full list of relevant references, has already 
appeared [^], we use this opportunity to give a less technical and more intuitive presentation. 

2. A brief review of the EK and QEK models 

The U (N) EK model of Ref. [||] is a matrix model whose connected correlation functions are 
expected to have, under certain assumptions, the same large-A^ limits as appropriately chosen cor- 
relation functions in a U (N) pure gauge theory. It is a specific example of an orbifold projection 
mapping a "mother" theory (the gauge theory in our case) to a "daughter" theory (here the matrix 
model), that, under certain conditions, becomes an equivalence in the large-A^ limit. This equiva- 
lence applies only to "neutral sectors" of the two theories, which consist of correlation functions 
invariant under translations in the gauge theory, and invariant under the U(\) d "center" symmetry 



in the matrix model [ |10| ] (d is the number of space-time dimensions). Obtaining the reduced model 
from the gauge theory is easy: one simply sets to zero all Fourier components of the gauge fields 
except for the zero mode and so performs the mapping U^(x) — > Up. Substituting this into the 
pure gauge Wilson action gives the partition function of the EK model, a single-site lattice gauge 
theory: 1 

Z EK = [dU exp{-S EK ) , S EK = -Nb £ 2ReTr (fytfytf^) . (2.1) 

The essence of large-Af reduction is that vacuum expectation values (VEVs) of neutral opera- 
tors that map into each other have coinciding large-A^ limits. For this to hold the vacuum has to be 
invariant under the symmetries defining the projection, for otherwise the desired expectation values 
do not lie in the neutral sectors. Indeed the original EK paper [||] focused on center-invariant VEVs 
of Wilson loops projected to zero momentum. Assuming implicitly that translation invariance is 
unbroken in the gauge theory, while stressing explicitly that the center symmetry must remain in- 
tact in both the gauge theory and matrix model, the authors proved that these VEVs become equal 
to corresponding quantities in the matrix model at large-A^. (Technically, they showed that the 
Dyson-Schwinger equations obeyed by Wilson loops in the two theories coincide.) 

Unfortunately, the vacuum of the EK model was shown to spontaneously break its center sym- 
metry for weak enough bare lattice coupling if d > 2 111], Oj. Consequently, the EK model 

'The U(l) d center symmetry is l/ M -> e' e "(/ M . The gauge symmetry is t/ M -> G Up G 1 ' for all /l,Ce U(N). 
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does not reproduce the neutral sector of large-A^ QCD. Initially this was seen in the weak-coupling 
limit by analyzing the effective potential, V e s, felt by the eigenvalues of the link matrices [@]. 
The eigenvalues are attracted, leading to a ground state with proportional to the identity up to a 
phase, which spontaneously breaks the center symmetry To overcome this they suggested quench- 
ing the eigenvalues, thereby defining a new prescription for calculating neutral-sector expectation 
values, ((^quenched- This prescription has two parts: 

1. Calculate \ ff) p for a frozen set of eigenvalues: Eig(f/^) = e iPa ; a = 1,2, . . . ,N. Here (,) p 
denotes an average with an action obtained from Eq. ( |2T] ) by freezing the fluctuations in 
the eigenvalues of U^. This means writing = e ip>1 Vn, with = diag(p^, . . . ,p^) the 
frozen eigenvalues and G U (N), and replacing DU^ with DV^. 

2. Average over the p" with a measure d n(p) that is both U(l) d invariant and distributes each 
of the p^ uniformly and independently over [0,2n) when N — > °°: 

) quenched = / Ua,^ d^(pa ) {&)p. 



We stress that this prescription differs substantially from the ordinary way one calculates expecta- 
tion values in field theory. In particular, for the quenched prescription to coincide with an ordinary 
field theory average one would need to use dpL(p) °< YlaudPa^ ~ v * B u') [13]. This is a highly non- 



uniform measure over the eigenvalues, as noted above. By removing the V e n factor, e.g. by using 
the uniform measure djx{p) oc JT a dp„, one expects the center symmetry to remain unbroken. 
It is illuminating to understand the quenching prescription in a different way 



approach is first formulated for two-index scalar fields. For example, one can define a mapping 
between the large- volume theory with N xN adjoint fields <p(x) and a theory of N x N matrices 
(j): (j) ab (x) — > e^^Pa-PbYH § ab . Here a,b are color indices, is the space-time coordinate, and 
Pa are predefined variables which are referred to as momenta for a reason that will become clear 
below. Substituting this mapping into the field theory action and observables yields the matrix 
model, whose expectation values depend on the p^ 2 Next, one is instructed to integrate such 
expectation values over p. Under certain assumptions, the result will become equal, when N — ► °°, 
to the corresponding expectation values in the field theory. 

This equality can be shown to any order in perturbation theory (perturbing in, say, the three- 
point coupling of the scalar theory) [|l5|]. At large N, planar diagrams dominate. For such diagrams, 
one finds that, prior to the p integral, the perturbative expression in the matrix model at a given 
order coincides with the momentum space integrand of the corresponding Feynman diagram in the 
field theory. In this correspondence the momentum of the field <p ab (q) is given by the difference 
Qab = Pa ~ Ph m tne matrix model. The p integration can then be identified as being over the 
momenta flowing in the planar diagrams, and one recovers the field theory integral. 

Thus the two steps in the quenching prescription can be thought of as (1) calculating the 
contribution to an expectation value (ff) p of a point p^ 1 at a Euclidean Brillouin Zone, and (2) 
integrating over the Brillouin zone uniformly to get the full (£?) q Uen ched- This way of embedding 
space-time (or rather its first Brillouin zone) into index space is typical of volume reduced models 
and was already noted in Refs. p], [T3|]. 

Performing the mapping in a lattice gauge theory, so that ty ab {x) becomes U^ b (x), one obtains 
a matrix model with an action similar to that in Eq. ([O]), but with replaced by e ,p>1 . The 
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momentum factor e' p can, however, be absorbed by a change of variables, so one ends up back at 



the problematic EK model. Various approaches have been used to avoid this problem. Refs. [ |14| ] 



restrict their study to 0(N) models, while Refs. [|T|, [16|] change the measure of the path integral: 



1. Ref. [0]: J U{N) DU^J U{N) DU^ J U{N) DV^8(U^ -V^e^V*) (up to /^-dependent terms), 
thus forcing the eigenvalues of U fl e ip to be e xp " . Intuitively, is being forced to fluctu- 
ate around unity, so that absorbing p 11 into £/„ is impossible. Note that the integral over Vn 
could equally well be over the coset U (N) /U(l) N in which multiplication from the right by 
aU(l) N matrix is divided out. 

2. Ref. [Q: Ju^DU^ — > Ju(n)/u(i) n ^h (together with a gauge-fixing term). Here one re- 
stricts the integration regime to the coset such that cannot absorb the e' p>1 factor. 

The first prescription is, in fact, identical to the QEK model of Ref. as can be seen by using 
the 5— functions to perform the U integrations [|l5|]. The second prescription is different and we 
call it the DW model. An obvious advantage of the QEK over the DW model is that the "reduced" 
U(N) gauge invariance U^e ipil — * QU ii e ip & is realized in the QEK model as V - ^ — > OV^. By 
contrast, the DW model has no gauge symmetry: if one writes QXJ^e ipll Q^ = U'^e 1 ^ then, in 
general, U'^^U (N) /U (1) N . Indeed, the DW model is defined including gauge-fixing so as to avoid 
this problem JT^]. This means that a numerical study of the DW model appears quite nontrivial. 

Let us restate in another way why one needs to change the measure of the path integral over the 
. In the weak coupling regime, b — > «>, configurations close to the minima of the action dominate. 
To avoid the symmetry-breaking of the EK model, one wants there to be a single minimum (up to 



gauge transformations). For the action ( |2.1[ ) with replaced with U jl e ip , however, there are 
multiple minima, occurring when the U^e ip ^ in all directions are arbitrary diagonal matrices (up to 
gauge transformations). By changing the measure as in Refs. [15, [16|] one aims to make the only 



minimum available that at = 1. This is indeed what happens, by construction, in the DW model. 
One might expect the same for the QEK model: the action is minimized when e ip ^ € U(\) N 
(up to gauge transformations), and this condition is satisfied by = 1. This, however, is not the 
whole story: as has long been known, there are other solutions to this equation given by = 
with a U (N) /U (l) N matrix that, when conjugating e' Ptt , gives rise to a permutation of the color 
indices of p" . These minima are nonperturbative (P^ is very far from the unit matrix) and we find 
that they play a central role in the (in)validity of large-Af quenched reduction. 

We began this section by describing reduction as an example of the orbifolding paradigm. 
Once one quenches the daughter theory, however, this paradigm does not apply in a straightforward 
way. This is because VEV's in the daughter QEK model are calculated in a way that separates the 
gauge fields into quenched and unquenched degrees of freedom — a separation that is absent from 
the mother gauge theory. Nonetheless, an argument for quenched reduction has been given using 



the Dyson-Schwinger equations for Wilson loops []15j]. These "loop equations" are different in 
the gauge theory and QEK model, but the differences are proportional to quenched expectation 
values of center-invariant quantities composed of center non-invariant Wilson loops. An example 
of such a quantity is (|Wc| 2 ) qU enched> where Wc is center non-invariant. Reduction can hold only if 



2 Up to right-multiplications by U(l) N matrices which we can avoid by restricting to U (N) /U(l) 
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such quantities vanish as N — > °°. They will vanish if quenched expectation values factorize, for 
then ( | Wc | 2 )quenched = I (Wc)quenched | 2 + 0(\/N), and the enforcement of center-symmetry in the 
QEK implies that \(Wc) quenched | = 0. The crucial question is thus whether factorization holds for 
quenched expectation values. 

3. (In) Validity of large-N quenched reduction 

In this section we discuss analytic, weak-coupling considerations that extend those presented 



in Refs.[]4|] and [ |1 lfl to include the effects of the multiple classical minima which, as discussed 
above, come in the form of permutations of the pa- To be precise, the minima are £/™"V pP = 
ymm g i>" y^ mm _ gtp" > w j t j 1 _ an( j <j( fl ) one f th e jyi permutations of the indices a. Here 
we describe how this classical degeneracy is removed by quantum fluctuations, with the resulting 
ground state corresponding to a certain permutation which is unrelated to the original choice of 
Pa - This, by itself, is a strong indication for the breakdown of quenched large-N reduction, since 
it invalidates the idea that for each value of p^, the value of (&) p is the contribution of a point 
corresponding to p^ in the Brillouin zone. 

To show how quantum fluctuations remove the classical degeneracy we minimize the effec- 
tive potential V e ff over the space of all permutations of the color indices of the input momenta pa 
(with permutations in different directions being independent). For brevity, we consider here the 
case where the ,/V momenta in each direction, pi} , are always a permutation of the "clock" mo- 
menta p c a {ock = 2na/N. The situation for other distributions of momenta is similar, as discussed in 
Refs. [0J. For clock momenta we find that, for d > 3, the minimum of V e ff occurs for permutations 
in which the momenta in all directions are ordered similarly. In equations, the statement is that 
p% —p v a = for all a, and we refer to this as "locking" of the momenta. There is thus an overall 
shift between different Euclidean directions, whose value is an arbitrary clock momentum that can 
depend on /i and v, but which is independent of the color index. This implies that the difference 
q\ = pa — p^, which we recall plays the role of the gluon momentum in Feynman diagrams, is 
independent of fl. This in turn means that instead of integrating over the full Brillouin zone (as 
would be the case if the p^ in different directions were independent) one is effectively integrating 
along the diagonal of the Brillouin zone, i.e. over a one-dimensional line! 

It is useful to have order parameters that can detect such "locking", and a simple choice is 
the center-non-invariant reduced Wilson loops M^ v = hxxU^Uv and _ v = 4tr U^ul? In- 
deed, when the locked ordering takes place one finds that _ v = e' aMV , while M^ iV = 0. By 
contrast, if no ordering takes place and all pa are random permutations of p" lock , then M^ ± v ~ 
0(\/N). In Fig. H we show the dependence, in d = 4, of Veff on the combined order parameter 
M = ^Lji (iM^vl + |-M/i,-v|)> where we have used a Monte-Carlo (MC) to sample permutations. 
One sees that Veff has a minimum at the "fully locked" ground state with M = 6, with states having 
little locking (M 1) having a free energy of 0(N 2 ) higher. 

We now come to a crucial observation: if locked ordering occurs then large-A^ factorization 
breaks down in the quenched theory. To see this, recall that the quenched prescription instructs us 
to integrate over the pH in a center-symmetry invariant way. Performing a center transformation on 
the pa (pa — > Pa + will lead to a center-transformed locked vacuum, since the action is center 



We note that with clock momenta one always has tr = 0. 
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Figure 1 : Right panel : The dependence of V t g /N 2 on the combined order parameter M for random permuta- 
tions of the clock momenta and for various values of N. Left panel: Scatter plot of the values of three chosen 
M/xv's in a run for N — 4-0 and 'clock' momenta. We refer to |Q| for technical details. 



invariant. In this transformed vacuum, a^ v — > a^ 11 + 8^ — 8 V , so the phase of the Mn V is shifted. 
Integrating over the 6^ , which is part of the integration over the p% , will thus lead to the vanishing 
of the order parameters, e.g. in the example above, (M /J) _ v )q Ue nched — J n da^ v e' a ^ v = 0. By 
contrast, if one calculates (\Mn ,-v| 2 )quenched> the phase a^ v cancels, and so _ v | 2 ) quenched ~ 
0(1) ^ | (M^ v ) quenc hed| 2 - This invalidates large-Af factorization in the quenched theory, which, as 
noted above, implies that the gauge theory and QEK model have different loop equations, and so 
reduction fails. 

The argument for locked ordering described above was perturbative. To check whether it 
occurs beyond perturbation theory we performed a detailed numerical study of the QEK model at 
intermediate and strong couplings using MC techniques for 20 <N < 200. We do not discuss any 
particulars of these calculations here, referring instead to Refs. [§] (the second reference includes 
an appendix describing the update algorithms we used, which are nontrivial since the model is 
quartic in the updated fields). 

In our numerical studies we measure the order parameters M^ v in the intermediate coupling 
regime, and find clear evidence for correlations in the eigenvalue orderings. For example, in the 
right panel of Fig. |] we show a typical scatter plot of data obtained for Mj.\ , M4.3, and M3 _2- This 
combines results from 20 independent input choices of the p^, , all being random permutations of 
^ciock Wri at we see from the figure is that each of these 20 MC simulations settled into a different 
vacuum characterized by a distinct value for the phase a^ v of M^ v . This is precisely what one 
expects when locking of the eigenvalue ordering occurs. 

We have also obtained direct evidence for the breakdown of large-A^ reduction: there are large 
discrepancies between the plaquette values and the structure of the phase diagram of the QEK 
model and the large-A^ gauge theory We have checked that these conclusions are insensitive 
to the precise form of the quenched eigenvalue distribution, d}i{p), and to the way we perform 
the quenched average. We also considered values of ,/V up to 200 to look for a late onset of 1 /N 
behavior, but find none. Therefore we conclude that the momentum quenched large-A^ reduction of 
SU (N) lattice gauge theories fails in the weak-coupling regime and thus in the continuum limit. 
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4. Possible future directions 

Our result, together with those of Refs. [B], imply that the only known single-site models that 
can reproduce the properties of QCD at large-TY are the following: (I) The "deformed" Eguchi 



Kawai of Ref. [17]; (II) the single-site obtained by adding adjoint fermions to the Eguchi-Kawai 



model [jg]; and (III) the momentum quenched model of DW Q16Q. For a brief discussion on the 
first two we refer to Refs. As far as we know, the third alternative has never been explored 
nonperturbatively; its obvious advantage over the standard type of momentum quenching is that its 
path integral has a unique minimum in weak coupling, with no difficulties related to permutations. 
As mentioned above, however, this model involves gauge fixing and so it is clearly more involved 
to study numerically. We leave the exploration of all these single-site models to future studies. 
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